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Abstract 

We discuss production of two pairs of cc in proton-proton collisions at the LHC. Both double- 
parton scattering (DPS) and single-parton scattering (SPS) contributions are included in the anal- 
ysis. Each step of DPS is calculated within A:j-factorization approach, i.e. effectively including 
next-to-leading order corrections. The conditions how to identify the DPS contribution are pre- 
sented. The discussed mechanism unavoidably leads to the production of pairs of mesons: DiDj 
(each containing c quarks) or DiDj (each containing c antiquarks). We calculate corresponding 
production rates for different combinations of charmed mesons as well as some differential distri- 
bution for [D^D^ + D^D^) production. Within large theoretical uncertainties the predicted DPS 
cross section is fairly similar to the cross section measured recently by the LHCb collaboration. 
The best description is obtained with the Kimber-Martin-Ryskin (KMR) unintegrated gluon dis- 
tribution, which very well simulates higher-order corrections. The contribution of SPS, calculated 
in the high-energy approximation, turned out to be rather small. Finally, we emphasize significant 
contribution of DPS mechanism to inclusive charmed meson spectra measured recently by ALICE, 
ATLAS and LHCb. 
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I. INTRODUCTION 



There has been recently renewed interest in studying double-parton scattering (DPS) 
effects in different reactions (see e.g. [1] and references therein). Very recently we have 
shown that the production of cccc is a very good place to study DPS effects [2]. Here, the 
quark mass is small enough to assure that the cross section for DPS is very large, and large 
enough that each of the scatterings can be treated within pQCD. The calculation performed 
in Ref. [2] were done in the leading-order (LO) coUinear approximation. This may not be 
sufficient when comparing the results of the calculation with real experimental data. In the 
meantime the LHCb collaboration presented new interesting data for simultaneous produc- 
tion of two charmed mesons [3] . They have observed large percentage of the events with two 
mesons, both containing c quark, with respect to the typical production of the corresponding 
meson/antimeson pair (o'DiDj /o'D.Dj ~ 10%), despite of the very limited LHCb acceptance. 

Is the large effect a footprint of double parton scattering? We wish to address the issue 
in this paper. In addition, we shall estimate cccc production via single-parton scattering 
(SPS) within a high-energy approximation [4]. This approach should be an efficient tool 
especially when the distance in rapidity between cc or/and cc is large. 

Another evidence for the DPS effects can be a missing cross section in the inclusive 
charmed meson distributions observed recently in Ref. [5]. The measured inclusive cross 
sections include events where two D (or two D) mesons are produced, therefore correspond- 
ing theoretical predictions should also be corrected for the DPS effects. 

In Ref. [6] the authors estimated DPS contribution based on the experimental inclusive 
D meson spectra measured at LHC which, as discussed in our paper, may be too crude 
approximation. In addition, in their approach fragmentation was included only in terms 
of the branching fractions for the transition c ^ D. In our approach we shall include full 
kinematics of hadronization process. Here we wish to show also first differential distributions 
on the hadron level to be confronted with recent LHCb experimental data [3] . 

II. THEORETICAL FRAMEWORK 

In the present analysis, when considering pp — )■ ccccX reaction, we concentrate primarily 
on double-parton scattering effects. In Section III.B we will show that the single-scattering 
contribution to double-charm production is much smaller, especially in the LHCb kinematics. 

A. Double-parton scattering 

In LO coUinear approximation the differential distributions for cc production depend 
e.g. on rapidity of quark, rapidity of antiquark and transverse momentum of one of them 
(they are identical) [2]. In the next-to-leading order (NLO) coUinear approach or in the 
/crfactorization approach the situation is more complicated as there are more kinematical 
variables necessary to describe the kinematical situation. In the /Cffactorization approach 
the differential cross section for DPS production of cccc system, assuming factorization of 
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the DPS model, can be written as: 
^^DPS^pp _^ ccccX) 



1 da^^^ipp — > ccXi) da^^^{pp ccX: 



^\ (2.1) 



2(7e// dyidy2d'^pi,td'^p2,t dy^dyid'^p^^td'^p^, 
When integrating over kinematical variables one obtains 

a^^^pp ccccX) = -^a^^^{pp ccX^) ■ a^^^{pp 00X2). (2.2) 

These formulae assume that the two parton subprocesses are not correlated one with each 
other. The parameter cTe// in the denominator of above formulae can be defined as: 

-1 



^eff 



d%{T{h)f 



(2.3) 



where the overlap function 



T{h) = j f{b,)f{h - b)d%, (2.4) 

if the impact-parameter dependent double-parton distributions (dPDFs) are written in the 
following factorized approximation [7, 8]: 

Tij{xi,X2;bi,b2;iJ,l,i^l) = Fij{xi,X2; i^j, iJ,l)f{bi)f{b2). (2.5) 

Experimental data from Tevatron [9] provide an estimate of (Tg// in the denominator of 
formula (2.2). Corresponding evaluations from the LHC are expected soon. In our analysis 
we take (Je// = 15 mb. In the most general case one may expect some violation of this simple 
factorized Ansatz given by Eq. 2.2 [8]. 

In our present analysis cross section for each step is calculated in the /cj-factorization 
approach, that is: 

dgSPSjpp ^ ccX^) 1 r d^htd^k2t 

dyidy2d'^Pi^td'^P2,t IGtt^P 7 tt tt 9 9 ^cc 

X 5^ (kit + ht - Pit - P2t) J^{xi, kf^, i/)J'{x2, kit, A«^)> 
da^P^{pp CCX2) 1 f d^k^td^kit- 



dysdy4d^P3^td^P4,t 16n^P J tt n '"'"'^ ^ 

X (kst + ku - P3t - P4t) J^{x3, kit, I^^)^{xa, kit, A*^)- (2-6) 

The matrix elements for g*g* cc (off-shell gluons) must be calculated including transverse 
momenta of initial gluons as it was done first in [10-12]. The unintegrated (fcrdependent) 
gluon distributions (UGDFs) in the proton are taken from the literature [13-15]. Due to 
the emision of soft gluons encoded in these objects, it is belived that a major part of NLO 
corrections is effectively included. This is in analogy to initial state parton shower in Monte 
Carlo generators and strongly depends on technical construction of UGDF (see Ref. [5]). The 
framework of the /ct-factorization approach is often used with success in describing inclusive 
spectra of D or B mesons as well as for theoretical predictions for so-called nonphotonic 
leptons, products of semileptonic decays of charm and bottom mesons [16-22]. 
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B. Single-parton scattering 

The total cross section for the production of cccc final state via single gluon-gluon inter- 
action can be calculated in the parton model approach as: 

a{pp — >■ CCCC] W"^) — j dxidx-2 g{xi, g{x2, //p) cr{gg — >■ cccc; XiX2W'^) . (2-7) 

Here g{x,fi^) is integrated (collincar) gluon distribution in a proton (PDF), and W is the 
proton-proton center of mass energy. In practice the integration is done in log^^o^i 
log^o^a, including the corresponding jacobian of transformation. The elementary cross sec- 
tion of Eq. (2.7) enters at s = XiX2W'^ > 16m^. The parton level cross section in (2.7) is 
therefore very useful in order to obtain differential distributions in invariant mass of the cccc 
system. 

In the present calculation we concentrate on LHC energies and consider the gg cccc 
subprocesses only. In the high-energy approximation the elementary cross section can be 
written in the compact form (see Ref. [4]): 

N'^ — 1 An^a"^ -> -> d'^ki d'^k2 d^q 

da{gg ^ cccc) = J^TJ^ h-^cc{zu h, W-(^2, h, -q) dzij^ "^^^J^J^- 

(2.8) 

Here the Ig^cc{zi,ki,qi) and Ig^cc{z2,k2,q2) factors, called impact factors, describe the 
coupling of pairs of cc associated with the first and second gluon/proton, respectively. Above 
zi and Z2 are longitudinal momentum fractions of quarks with respect to parent gluons in 

— * 

the first and second pair, respectively, and ki their respective transverse momenta, q is 
exchanged transverse momentum and hq is gluon mass which can be put to zero at least 
mathematically. At low energies this formula must be corrected for threshold effects [4]. The 
differential cross sections for pp ccccX can be obtained by replacing the a{gg cccc) by 
da{gg — )■ cccc) in Eq.(2.7). Details about how the arguments of ag are chosen are discussed 
in Ref. [4]. 

Our approach includes subprocesses coherently to be contrasted to Ref. [6] where they 
were separated one from each other to simplify calculations. In addition, we get a practical 
agreement with results of calculations in Ref. [23]. 



C. Double meson production 

Kinematical correlations between quarks and antiquarks are not accessible experimen- 
tally. Instead one can measure correlations between heavy mesons or nonphotonic electrons. 
In this paper we will analyze kinematical correlations between charmed mesons. In particu- 
lar, we are interested in correlations between Di and Dj mesons (both containing c quark) 
or between and Dj mesons (both containing c antiquark). In order to calculate correla- 
tions between mesons we follow here the fragmentation function technique for hadronization 
process: 

da{pp DDX) _ f D^_^d{zi) _ Dc^d{z2) _ da{pp ccX) 
dyidy2d^p^^d'^p§^ ^ J Zi Z2 dyidy2d^p'i td^p^^^ ^ ^ ' 
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where: pl ^ = pg,* = ^ and meson longitudinal fractions zi, Z2 G (0, 1). We have made 
approximation assuming that yi, y2 and (f) are unchanged in the fragmentation process. The 
multidimensional distribution for both c quarks (or both c antiquarks) is convoluted with 
fragmentation functions simultaneously for each of the two quarks (or each of the two anti- 
quarks) . As a result of the hadronization one obtains corresponding two-meson multidimen- 
sional distribution. In the last step experimental kinematical cuts on the distributions can 
be imposed. Then the resulting distributions can be compared with experimental ones. For 
numerical calculations here we apply often used in the case of heavy quarks, the Peterson 
fragmentation function [24]. We have shown in Ref. [5] that this scheme works very well in 
the case of inclusive meson spectra as well as for D^D'^ kinematical correlations. 

III. RESULTS 
A. Parton level 

We start from inclusive distributions of charm quarks (or antiquarks). As discussed in 
Ref. [5] the standard single-parton scattering contribution to pp — > ccX seems insufficient 
to describe inclusive spectra of charmed mesons as measured by the ATLAS, ALICE and 
LHCb collaborations [25-27]. The cccc production also contributes to the inclusive charm 
production. In Fig. 1 we show such a contribution to transverse momentum distribution 
(left panel) and rapidity distribution (right panel) together with theoretical uncertainty 
band. In this calculation the Kimbcr-Martin-Ryskin (KMR) UGDF [13] was used with 
the MSTW08 [28] collinear gluon PDF. The solid line corresponds to central value of our 
predictions. The uncertainties are obtained by changing charm quark mass rric = 1.5 ± 0.3 
GeV which in general is not well known and by varying renormalization and factorization 
scales [i^ — iiji — fjLp — (mf, where ( e (0.5; 2). The shaded bands represent these both 
sources of uncertainties summed in quadrature. As a reference point we plot contribution 
from standard single-scattering cc production, obtained in the /c^-factorization approach 
(long-dashed line) as well as calculated with the help of FONLL code [29] (dash-dotted 
line). As can be seen both of these models are consistent and give very similar numerical 
results. It suggests that in the case of charm quark production the /cffactorization approach 
with the KMR UGDFs very well reproduces NLO corrections. These aspects of cc production 
were discussed in more detail in Ref. [5]. 

Since the DPS uncertainty band is very broad it becomes clear that this contribution 
is quite sizeable and must be included in the total balance of charm quark (atniquark) 
production. For comparison we show also DPS result obtained previously in Ref. [2] in the 
LO collinear approach. It is much smaller than the /cj-factorization result, especially at 
larger transverse momenta. 

In Fig. 2 we compare DPS results for transverse momentum (left panel) and rapidity 
(right panel) distributions obtained with different UGDFs from the literature [13-15]. The 
KMR UGDF gives the largest cross section. Numerical results of DPS are more sensitive to 
the choice of UGDFs than in the case of SPS cc production, which can be understood by 
different power of UGDFs in the cross section formula (fourth in DPS cccc versus second in 
SPS cc). We use here also the KMS [14] and Jung setA+ [15] paramctrizations. In turn, in 
Fig. 3 we confront theoretical uncertainties of SPS single pair (cc) and DPS two-pair {cccc) 
production. Again uncertainty of the two pair production is much larger than that for single 
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FIG. 1: Transverse momentum (left) and rapidity (right) of charm quarks from SPS cc (long- 
dashed line) and DPS cccc (solid line with shaded band) production. In this calculation the KMR 
UGDF was used and the factorization scale and quark mass for the DPS contribution were varied 
as explained in the figure. For comparison LO collinear DPS distribution (dotted line) and FONLL 
SPS cc result (dash-dotted line) are shown. 
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FIG. 2: Transverse momentum (left) and rapidity (right) distributions of charm quarks produced 
in DPS cccc production for different unintegrated gluon distributions. 



In Ref. [2] we have proposed several correlation distributions to be studied in order to 
identify the DPS effects. Here we present the same distributions as in Ref. [2] but within 
the fc^-factorization approach. In Fig. 4 we show distributions in invariant mass M^s (left 
panel) and rapidity difference of quarks/ ant iquarks Idj// = Vc — Vc (right panel) from the 
same scattering (ciC2 or C3C4) and from different scatterings (ciQ or C3C2 or C1C3 or C2C4) for 
various UGDFs specified in the figure. The shapes of distributions in the figure are almost 
identical as their counterparts obtained in LO collinear approach in Ref. [2]. 

In Fig. 5 we show distributions in azimuthal angle difference between quarks/antiquarks 
(fee from the same and from different scatterings. While in the case of the same scattering 
distribution strongly depends on the choice of UGDF the quarks/antiquarks from different 
scattering are not correlated which is inherent property of the simple factorized model. Our 
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FIG. 3: Comparison of the SPS cc and DPS cccc contributions to the inclusive charm quark 
(antiquark) production together with theoretical uncertainties due to the choice of scales and those 
related with quark mass (summed in quadrature). 




FIG. 4: Distribution in the invariant mass of quark/atiquark Mcc (left) and distribution in the 
rapidity distance between quarks/antiquarks Idi// (right) from the same (ciC2 or C3C4) and from 
different scatterings (ciq or C3C2 or C1C3 or C2C4), calculated with different UGDFs. 



distinguishing of scatterings can be done only in the model calculation. Experimentally 
one observes both types together after hadronization which naturally may bring additional 
decorrelation. 

Finally, we present distribution in transverse momentum of the pair of quarks pf. In 
LO coUinear approach the distribution for emission in the same scattering is very different 
from the case of emissions from different scatterings [2]. The picture in the fcrfactorization 
approach is, however, very different. The respective distributions for the same and different 
scatterings are rather similar. This means that transverse momentum of the pair may not 
be the best quantity to be used in order to identify the DPS effects. 
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FIG. 5: Distribution in azimuthal angle (/Jcc between quarks/antiquarks from the same scattering 
(left) and from different scatterings (right), calculated with different UGDFs. 
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FIG. 6: Distribution in transverse momentum of the quark/antiquark pair from the same 
scattering (left) and from different scatterings (right), calculated with different UGDFs. 



B. Meson level 

Production of two pairs of cc on the partonic level leads to the situations that very often 
two mesons, both containing c quarks or/and both containing c antiquarks, are produced 
on the hadronic level in one event. Therefore the presence of two such mesons may be 
considered as a signal of production of cccc on the partonic level. Recently, the LHCb 
collaboration performed a first measurement of D^Dj + DiDj production in the fiducial 
range of the detector acceptance 2 < ya < 4 and 3 < < 12 GeV [3]. As described in 
Section II we have prepared a code which keeps track of full kinematical information about 
each of both quarks or both antiquarks and similar information about both mesons. Such a 
multidimensional map is used then to impose adequate experimental cuts. 

In Table 1 we have collected DPS cross sections for different pairs of mesons relevant 
for considered kinematics obtained with different unintegrated gluon distributions. As was 
shown in Ref. [5] theoretical predictions for production of charmed meson pairs in the LHCb 
kinematics are very sensitive to the value of Ec parameter in the Peterson fragmentation 
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TABLE I: Total cross sections for meson-meson pair production for three different UGDFs. 



Mode (rgfP [nh] KMR +(mc) Jung setA+ KMS 

ec = 0.05 ec = 0.02 ec = 0.05 ec = 0.02 = 0.05 £c = 0.02 
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function. There, rather harder functions (with smaller Ec) are suggested for better descrip- 
tion of experimental data, which is also in agreement with observations made in the FONLL 
framework [30, 31]. Therefore we present results for two different values of the Ec parameter. 
Here we have added together cross sections for charge conjugated channels: aOiDj + <^DiDj- 
The calculated cross sections are somewhat smaller than the experimental ones. Only the 
upper limit of our predictions with the Kimber-Martin-Ryskin UGDF and with Ec = 0.02 in 
the Peterson fragmentation function gives results which are close to the experimental data, 
taken the uncertainties on the choice of the factorization/renormalization scale and on the 
charm quark mass. 

So far we have considered only DPS contribution to DiDj (or DiDj) production. In 
Fig. 7 wc show in addition corresponding SPS contribution. The SPS contribution to the 
transverse momentum distribution (left panel) is more than two orders of magnitude smaller 
than the DPS one. For the rapidity distribution (right panel) the difference is only one order 
of magnitude. This effect is shghtly unintuitive. However, it can be understood by a com- 
parison of the two-dimensional distributions in rapidity of one and the second D meson for 
DPS and SPS production (see Fig. 8). In the case of DPS the two mesons are not correlated 
(in this plane), in contrast to the SPS mechanism, where they are strongly anticorrelated. 
When one meson is produced in forward rapidity region the second is preferentially produced 
in backward rapidity region, or vice versa. One can also conclude that in the case of DiDj 
(DiDj) pair production, the specific LHCb kincmatical range leads to a dumping of the SPS 
cross section. The requirement that both D mesons have to reach the detector makes the 
SPS contribution almost negligible. Quite different conclusions can be drawn in the case of 
inclusive D meson measurements. 

In the present paper we have calculated SPS cccc contribution in the high-energy ap- 
proximation which may not be the best approximation for the LHCb kinematics where the 
distance between both c or both c is rather small. Therefore, to drawn definite conclusions, 
future studies of the pp — > ccccX process are needed and they must include a complete 
set of diagrams for the SPS cccc mechanism. Furthermore, if the improved calculations 
of SPS mechanism will not provide somewhat better description of the total cross sections 
measured by LHCb, one has to look for other mechanisms which can contribute and fill 
predicted missing strength. 

The LHCb collaboration presented also several differential distributions for the simulta- 
neous production of two DD and DD mesons. Here we consider only examples for D^D^ 
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FIG. 7: Distributions in transverse momentum (left panel) and rapidity (right panel) of single 
meson from the D^D^ pair events. The solid lines corresspond to DPS mechanism and the 
long-dashed lines represent contributions from SPS production of D^D^ pairs. Here we impose 
kinematical cuts adequate for the LHCb kinematics. 




FIG. 8: Two-dimensional distributions in rapidity of one meson and rapidity of the second 
meson for DPS (left) and SPS (right). 

(identical to D^D^) channel. 

In Fig. 9 we present distribution in transverse momentum of one of the mesons, 
provided that both are measured within the LHCb experiment coverage specified in the figure 
caption. Our theoretical distributions have shapes in rough agreement with the experimental 
data. The shapes of the distributions are almost identical for different UGDFs used in the 
calculations (left panel) and are independent on the choice of scales in the case of the KMR 
model (right panel). 

In Fig. 10 we show distribution in the D^D^ invariant mass M^oj^o for both D^^s mea- 
sured in the kinematical region covered by the LHCb experiment. Here the shapes of the 
distributions have the same behavior for various UGDFs and are insensitive to changes of 
scales as in the previous figure. The characteristic minimum at small invariant masses is a 
consequnce of experimental cuts (see Ref. [5]) and is rather well reproduced. 

Finally in Fig. 11 we show distribution in azimuthal angle fooDO between both D°'s. 
While the theoretical DPS contribution is independent of the relative azimuthal angle, there 
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FIG. 9: Transverse momentum distribution of D mesons from the D D pair contained in the 
LHCb kinematical region. The left panel shows dependence on UGDFs, while the right panel 
illustrates dependence of the result for the KMR UGDF on the factorization/renormalization scales. 




FIG. 10: M£)Q£)0 invariant mass distribution for D^D^ contained in the LHCb kinematical region. 
The left panel shows dependence on UGDFs, while the right panel illustrates dependence of the 
result for the KMR UGDF on the factorization/renormalization scales. 



is some small residual dependence on azimuthal angle in experimental distribution. This 
may show that there is some missing mechanism which gives contributions both at small and 
large A(f. However, this discrepancy may be also an inherent property of the DPS factorized 
model which does not allow for any azimuthal correlations between particles produced in 
different hard scatterings. We wish to emphasize in this context that the angular azimuthal 
correlation pattern for D^D^, discussed in Ref. [5], and for D^D^ (D^D^), discussed here, are 
quite different. The distribution for D^D^ {D^^D^) is much more flat compared to the D^D^ 
one which shows a pronounced maximum at (Puojyo = 180° (mostly from pair creation) and 
fjjojyo = 0° (mostly from gluon splitting) [5]. This qualitative difference is in our opinion 
a model independent proof of the dominance of DPS effects in the production of D^D^ 
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FIG. 11: Distribution in azimuthal angle (Pdo^o between both D^'s. The left panel shows depen- 
dence on UGDFs, while the right panel illustrates dependence of the result for the KMR UGDF 
on the factorization/renormalization scales. 



C. DPS cccc production and inclusive charmed meson distributions 

Since the DPS cross section is very large it is also very important to look at the DPS cccc 
contribution to inclusive charmed meson spectra. Let us consider for example transverse 
momentum distribution of a charmed Di meson. The corresponding DPS cccc contribution 
can be written as: 



dpt dpi,f 

, da^ , 

+ PA: 1 - ^A: h P2,t=P* < < 00, -2.1 < 7^2 < 2.1 

da^ 

+ Pd,Pd,-\ |pi,=pi(-2.1 <rii< 2.1,-00 < r/2 < 00) 

da^ , , N / X 

+ Pd^Pd-, — Ut=pt(-oo < ?7i < 00, -2.1 < ?72 < 2.1). (3.1) 



In the formula above Pd is a shorthand notation for the branching fraction Pc^d and 
(7^ is the cross section for D-mesons assuming artificially the branching fraction equal to 1. 
The formula above can be somewhat simplified when combining similar terms. 

In Fig. 12 we show inclusive one pair (long-dashed line), inclusive DPS two-pair con- 
tribution (dotted line) and the sum of both terms to transverse momentum distribution 
of different D mesons (solid line). The DPS cccc contribution is of the same order as the 
standard traditional SPS cc contribution. This is a completely new situation compared to 
what it was at smaller energies. The sum of both contributions almost describes the differ- 
ent experimental data. As discussed in the previous section the SPS cccc contribution can 
be of the order of 10% of the DPS cccc contribution. At higher energies one could expect 
even relatively larger DPS cccc contribution. A problem could start, however, that then one 
enters the region of really small gluon longitudinal momentum fractions x < 10~^ for which 
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P, (GeV) p, (GeV) 




FIG. 12: Inclusive transverse momentum distributions of different charmed mesons measured by 
different groups at the LHC. The long-dashed line corresponds to the standard SPS cc production 
and the dotted line represents the DPS cccc contribution. 



the gluon UGDFs (or PDFs) are not well known. In this case realistic models of UGDFs 
are badly needed. Do we have such a distribution at present? 
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IV. CONCLUSIONS 



In this paper we have discussed production of cccc in the double-parton scattering (DPS) 
and single-parton scattering (SPS) in the gg — )■ cccc subprocess. The double-parton scatter- 
ing is calculated in the factorized Ansatz with each step calculated in the /cj-factorization 
approach, i.e. including effectively higher-order QCD corrections. 

The cross section in the /cj-factorization approach turned out to be much larger than its 
counterpart calculated in the LO coUinear approach. The distribution in rapidity difference 
between quarks/antiquarks from the same and different scatterings turned out to have similar 
shape as in the LO coUinear approach. The same is true for invariant masses of pairs of 
quark-quark, antiquark-antiquark and quark- ant iquark, etc. The distribution in transverse 
momentum of the pair from the same scattering turned out to be similar to that for the 
pairs originating from different scatterings. 

We have calculated also cross sections for the production of Dj^Dj (both containing c 
quark) and DiDj (both containing c antiquark) pairs of mesons. The results of the calcula- 
tion have been compared to recent results of the LHCb collaboration. 

The total rates of the meson pair production depend on the unintegrated gluon distribu- 
tions. The best agreement with the LHCb result has been obtained for the Kimber-Martin- 
Ryskin UGDF. This approach as discussed already in the literature effectively includes 
higher-order QCD corrections. 

As an example we have also calculated several differential distributions for D^D^ pair pro- 
duction. Rather good agreement has been obtained for transverse momentum distribution 
of {D^) mesons and D^D^ invariant mass distribution. The distribution in azimuthal 
angle between both D°'s suggests that some contributions may be still missing. The single 
parton scattering contribution, calculated in the high energy approximation, turned out to 
be rather small. This should be checked in exact 2 — > 4 parton model calculations in the 
future. 

We have shown that the DPS mechanism of cccc production gives a new significant 
contribution to inclusive charmed meson spectra. For instance the description of the inclusive 
ATLAS, ALICE and LHCb data is very difficult in terms of the conventional SPS (cc) 
contribution [5]. 

Since we have shown that the DPS mechanism gives significant contribution to inclusive 
spectra of charmed mesons the estimate of DPS effects, presented in Rcf. [6] and based on 
experimental inclusive cross section, leads to an overestimation of the DPS effect. 

Summarizing, the present study of cccc reaction in the /ct-factorization approach has 
shown that this reaction is an extremely good testing ground of double-parton scattering 
effects. The LHCb kinematics is not the best in this respect. Both ATLAS and CMS 
collaborations could measure the production of pairs of DiDj and/or DiDj mesons with 
large rapidity distance where the DPS mechanism is predicted to clearly dominate over 
the SPS mechanism. Another potentially interesting place to investigate DPS effect is the 
pp — )■ J /ipJ /ipX reaction [32]. Similarly as for pp -> ccccX discussed here, the large rapidity 
gap between two J/V^'s should select clear sample of DPS mechanism. 
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